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High accuracy power spectra including baryonic physics in 
dynamical Dark Energy models 
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ABSTRACT 

The next generation mass probes will obtain information on non-linear power 
spectra P{k, z) and their evolution, allowing us to investigate the nature of Dark En- 
ergy. To exploit such data we need high precision simulations, extending at least up to 
scales of A: c± lO/i Mpc~^, where the effects of baryons can no longer be neglected. In 
this paper, we present a series of large scale hydrodynamical simulations for ACDM 
and dynamical Dark Energy (dDE) models, in which the equation of state parameter 
is z-dependcnt. The simulations include gas cooling, star formation and Supernovae 
feedback. They closely approximate the observed star formation rate and the observa- 
tionally derived star/Dark Matter mass ratio in collapsed systems. Baryon dynamics 
cause spectral shifts exceeding 1 % at A: > 2 — 3/iMpc~^ compared to pure n-body 
simulations in the ACDM simulations. This agrees with previous studies, although we 
find a smaller effect (^ 50%) on the power spectrum amplitude at higher fc's. dDE 
exhibits similar behaviour, even though the dDE simulations produce^ 20 % less stars 
than the analogous ACDM cosmologies. Finally, we show that the technique intro- 
duced in Casarini et al. to obtain spectra for any u;(z) cosmology from constant-w 
models at any redshift still holds when gas physics is taken into account. While this 
relieves the need to explore the entire functional space of dark energy state equations, 
we illustrate a severe risk that future data analysis could lead to misinterpretation of 
the DE state equation. 

Key words: galaxies: haloes - cosmology:theory, dark matter, gravitation - methods: 
numerical, N-body simulation. Hydro simulations 



1 INTRODUCTION 

There can be little doubt that Dark Energy (DE) is a nec- 
essary ingredient for any cosmological model approaching 
data (see, e.g., the data fits in WMAP5 release, Komatsu 
et al. 2009). In turn, the quest for the nature of DE is per- 
haps the main challenge in contemporary physics in spite of 
the fact that no available data confiict with ACDM models, 
in which the DE state parameter is ui = —1, as is true for 
a false vacuum. A false vacuum interpretation of DE, how- 
ever, leads to an unnatural value for its density (fine tuning 
problem) and to a number of coincidences characterizing our 
epoch and the primeval fiuctuation ampfitude (coincidence 
problem) . 

Many efl^orts have therefore been devoted to finding vi- 
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able alternatives. Among them, the possibility that DE is 
a scalar field, </>, self- interacting through a potential V[(j>) 
(Wetterich 1989, Ratra & Peebles 1989), has been widely 
explored. In this case, both the ratio between kinetic and 
potential energies, ij)^/2V{4>) and the DE state parameter 
w{z), depend on redshift, z. In this paradigm, observational 
data on 10(2) would yield V{(j}). We will refer to this class of 
models, as dynamical DE (dDE). 

Future large tomographic shear surveys, such as the Eu- 
clid project (see, e.g., Refregier et al. 2010), will measure 
matter density fluctuations and their evolution from the lin- 
ear to the non-linear regime with unprecedented precision, 
approaching ~ 1% (Huterer & Takada 2005). To translate 
these data into information on DE, we need to predict power 
spectra and the evolution of power spectra with similar pre- 
cision for a wide set of 'w{z) laws. The resulting constraints 
on w{z) will be a substantial step toward understanding DE 
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physics (see, e.g., Hu & Tegmark 1999, Dodelson & Zhang 
2005, Manera & Mota 2005, Mota 2008, La Vacca et al. 2008, 
2009, Kristiansen et al. 2010.). 

The evolution of the baryonic component of the Uni- 
verse is highly uncertain and can affect predictions for lens- 
ing observables, manifesting as modified structure growth 
for a fixed cosmic distance scale (Zentner et al 2008, Hearin 
& Zentner 2009). 

CoUisionless N-body simulations have been exten- 
sively used to obtain non-linear power spectra (e.g. Smith 
et al. 2003). These simulations have ignored the effect of 
baryons, which comprise ~ 15 % of the matter content of the 
observed Universe, (Komatsu et al. 2009). Although baryons 
make up a minor fraction of the Universal matter content 
and is distributed exactly the same as Cold Dark Matter 
(CDM) at the onset of the non-linear growth, the final dis- 
tribution of baryons in haloes is significantly different than 
CDM, because of its separate later dynamical evolution. 

Jing et al. (2006) have shown that, in the spectra of 
ACDM models, the difference between N-body simulations 
and simulations including baryon physics exceeds 1% at 
k > 2-3/iMpc~^, and increases to ~ 10-20% on smaller 
scales when k approaches lO/iMpc"^. Rudd et al. (2008) 
and Guillet et al. (2010) recently confirmed the need for 
simulations including baryonic physics, even though their 
results differed from Jing et al. (2006). 

The discrepant results originate from the different pre- 
scriptions used for star formation, SN feedback, AGN effects, 
and possibly, the numerical approach. To compare different 
wi^z) laws, these uncertainties need to be under control. Us- 
ing the same prescriptions for all models is essential to com- 
pare their spectra, as we do in this work. Their being tuned 
to observations is also a need, when aiming to fit future 
data. The prescription used here are tuned on available ob- 
servations, such tuning being however carried on within the 
frame of ACDM cosmologies. This is to be borne in mind, 
in view of model dependent effects that we shall detect, e.g. 
for what concerns star production in haloes. 

Dynamical Dark Energy (dDE) simulations, with a vari- 
able state parameter w[a) deduced from scalar field poten- 
tials admitting a tracker solution, have been performed since 
2003 (e.g. Klypin et al. 2003, Linder & Jenkins 2003, Dolag 
et al. 2004, Maccio et al. 2004, Solevi et al. 2005). These 
have been compared with constant-iu simulation outputs. 
Observables considered in these papers, only marginally in- 
cluded spectra. Recently Francis et al. (2007) have shown 
how spectral predictions for constant-u; models at z = 
can also be used to fit spectra of cosmologies with a state 
parameter given by a first degree polynomial: 



•][z) = Wo + w (1 — a) 



(1) 



More precisely, the spectrum P{k,0), in an assigned (A) 
model with state parameter 'w{z), can be approached by a 
suitable auxiliary constant-u; model (W) if: (i) in A and 
W, ^b,m,tot, h and as are equal, (ii) the constant DE state 
parameter w of W is tuned so that W and A have the same 
comoving distance from the Last Scattering Band (LSB) and 
2 = 0. In that case, spectral discrepancies stay < 1%, up 
to A; ~ 2-3/iMpc~^ (symbols here above have their usual 
meaning) . 

Francis et al. (2007) also tried extending their results 
to higher redshifts, but discrepancies between A and W in- 



creased to several percent; not enough to exploit forthcom- 
ing data. The required precision was, however, recovered 
through a new technique introduced by Casarini, Maccio 
& Bonometto (2009, Paper I hereafter) and tested with N- 
body simulations. 

The final aim of this work is to extend the techniques 
use in Paper I to scales that only hydrodynamical simulation 
can test. This is a serious challenge for the description of 
physical processes included in hydrodynamical algorithms. 
We will devote a part of our work to understanding these 
challenges. 

The focus of this paper is twofold: (i) Study the impor- 
tance of the effects of baryons on the matter power spectrum 
on the small scales that will probed by future surveys, bear- 
ing in mind that the physical assumptions will have to be 
carefully calibrated to approach the unprecedented preci- 
sion requirements, (ii) Extend the commonalities found in 
A''-body simulations in Paper I among spectra of different 
models down to scales where baryon physics cannot be dis- 
regarded. Accordingly, we shall use the same baryon physics 
in a series of simulations to compare the spectra of different 
dDE models. 

The paper is therefore organized as follows: In §2, we 
describe the approach used in Paper I. In §3, we motivate 
the choice of the models considered. §4 is devoted to de- 
scribing our simulations and the techniques used to analyse 
them. In §5, we present our hydrodynamical simulations. §6 
is subdivided in subsections that concern various aspects of 
the power spectra analysis. In §6 we test: (i) the relation 
between N-body and hydrodynamical spectra (ii) the effect 
of box size and numerical resolution (iii) the level at which 
model parameters can be discriminated; and, finally, (iv) the 
success of the extension of Paper I approach to the higher k 
range. A general discussion of the results closes our paper. 



2 APPROACHING DYNAMICAL DE 
SPECTRA 

Here is a short summary of the technique presented in Paper 
I, which improved the results of Francis et al. (2007) at z > 
0. Given an assigned model A, we introduce an auxiliary 
model W(z), for each z. A and W(z) are required to share 
the values of both ujb,c.m = ^b.c,mh^ and erg, at the selected 

2 . 

The first condition is satisfied at any redshift using the 
relationship between H and critical density pcr 



H"^ = (87rG/3)pcr 



(2) 



where H is the Hubble parameter. Multiplying both sides of 
this equation by Q.m (or fif,, Jlc) produces 



Q.mH 



(87rG/3)p„ 



(3) 



The r.h.s. of this equation scales as a~ , independent of 
model, since Um oc QmH^ (or ujb, uJc). Therefore, once the A 
and W models share ujb.cm at 2 = 0, they also share ijJb,c,m 
at every other 2. Accordingly, each yV(z) model shares the 
same u!b,c,m values. 

Conversely, the evolution of erg depends on the state 
equation of DE and its value at 2 = and can be obtained 
only from the constant DE state parameters 10(2), which is 
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Table 1. SIMULATION PARAMETERS 
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A.A Wg=-0.8 w = 
nW(z = 0.0) ACDM 
OW(z = 0.5) w=-1.035 
0'W{z=\.Q) w=- 1.079 
OMz = 2.0) w=-1.142 

WMAP5 

WMAP7 




-1.6 



■1.4 



■1.2 



-0.8 



-0.6 



W„ 



Figure 1. WMAP7 (solid curves) & WMAP5 (dashed curves) 
likelihood ellipses on the Wo-w' plane and the position of our 
models on that plane. The significant shift from WMAP5 to 
WMAP7 is not due to new CMB data, but mostly from the im- 
proved HST determination of Hq . Dotted lines yield the best like- 
lihood model. The selected model lies near the 1— cr contour for 
WMAP7, and within 2-cr's, for WMAP5 as well. Notice how much 
closer the W models (empty polygons) are than the assumed cos- 
mology (filled triangle). This highlights a serious observational 
risk (see JTJ. 



fixed by requiring that W{z) and A have equal distances 
between z and the LSB. 

The choice of Hq (the Hubble parameter at z = 0) is 
unconstrained, so we set it equal to -ffo in A to allow the 
simulated volumes to have the same side L in both Mpc and 
h~^Mpc units. Notice that a simple-minded generalization of 
Francis et al. (2007) criterion to high z would require equal 
iib,c,m{z) and, thence, H{z), and creates serious problems 
of sample variance and model comparison, at most redshift 
values. 

In Paper I, the mapping between yV{z) and A spec- 
tra was tested, using N-body simulations, up to fc ~ 2- 
3/iMpc~^. The discrepancies were mostly in the per- mil 
range. The techniques used in Paper I worked better for 
at higher z values. 



3 COSMOLOGICAL MODELS 

To extend our results to hydrodynamical simulations, we 
use three cosmological models, each of which is consistent 
with both WMAP5 and WMAP7 data. Two models share 
most cosmological parameters: baryon and CDM density 
fib — 0.046, Qc = 0.228, Hubble parameter (in units of 
100 km/s/Mpc) h = 0.7, r.m.s matter fluctuation ampli- 
tude at z; = on 8 /i~^Mpc scales, erg = 0.81, and primeval 
spectral index n — 0.96. The models differ in the assumed 
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dark energy equation of state, either w = ~1 for ACDM or 
using equation ((TJ with Wo = —0.8, w' — —0.754 for dDE 
(also sometimes referred to as the A model). We compare 
the dDE with various auxiliary constant-w models yV{z) . We 
plan to extend the tests performed on dDE to other 'w{z) 
expressions in future work, but there is currently no reason 
to believe that our dDE model is peculiar. In Figure [l] the 
A model and the auxiliary models are set on the WMAP 
likelihood ellipses on the Wo-w' plane. 

In order to extend the model comparison, we also con- 
sidered another ACDM model (ACDMI), with a greater 
CDM density Sic = 0.254, which yields a matter density 
parameter D,„, = 0.3. Thus, ACDMI is ~ 2 cr's from ACDM 
using the WMAP5 data. ACDM is run in boxes of differ- 
ent sizes to test the effectiveness of the mass and spatial 
resolution used. A summary of models and boxes used in 
this paper is presented in Table [T] ACDM coincides with 
W(z = 0). 



4 SIMULATIONS AND THEIR ANALYSIS 

We run both N-body and hydrodynamical simulations (here 
below also dubbed DMO and DMG, representing "DM- 
only" and "DM & gas" respectively). For DMO simulations 
we use PKDGRAV (Stadel 2001) that has been modified to 
enable variable w{z) as described in Paper I. For the hydro- 
dynamical simulations, we use GASOLINE, a multi-stepping, 
parallel Treesmoothed particle hydrodynamic (SPH) TV- 
body code (Wadsley et al. 2004). 

GASOLINE includes radiative and Compton cooling for a 
primordial mixture of hydrogen and helium. The star for- 
mation algorithm is based on a Jeans instability criteria 
(Katz 1992), but simplified so that gas particles satisfying 
constant density and temperature thresholds in convergent 
flows spawn star particles at a rate proportional to the local 
dynamical time (see Stinson et al. 2006). The star formation 
efficiency was set to 0.05 based on simulations of the Milky 
Way that satisfied the Kennicutt (1998) Schmidt Law. The 
code also includes supernova feedback as described by Stin- 
son et al. (2006), and a UV background following Haardt 
& Madau (1996); see Governato et al. (2007) for a more 
detailed description of the code. We applied the same mod- 
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Figure 2. Mass functions for the ACDM model at two dif- 
ferent redshifts {z = 0; 0.8) in the N-body and Hydrodynam- 
ical simulations: DMO and DMG respectively. Results for the 
L = 64/i~ Mpc and L = 256h~ Mpc boxes are shown with 
different symbols. Dashed and dotted lines show the Sheth & 
Tormen (2002) prediction for 2 = and z = 0.8. 



ifications to GASOLINE that we implemented in pkdgrav in 
Paper I to handle the dDE models. 

All of the simulations use 256^ CDM particles (and 
256'^ gas particles) and most of them are boxes 256 h~^Mpc 
on a side. For the ACDM, dDE models, and the related 
auxiliary models, the CDM particles each have a mass 
of mch/MQ = 7.61 x lO"' in the DMO simulations. In 
the DMG simulations, the CDM particles have a mass of 
irich/MQ = 6.33 x 10^" and the gas particles have a mass of 
rubh/MQ = 1.28 x 10^°. For the ACDMl model, rUch/MQ = 
8.33 X IQi" in the DMO and mch/MQ = 7.05 x 10^°, 
rubh/MQ = 1.28 x lO" in the DMG. 

The force resolution (softening) is 1/40 of the intra- 
particle separation. For L — 256/i'-iMpc, this corresponds 
to a distance e ~ 25 h^^kpc (wavenumber k. = 27r/e ~ 
150 h Mpc~i). This softening provides spectral resolution up 
to A: = lO/iMpc^^. Each simulation is started at Zin = 24 us- 
ing the same random realization of the density field. In order 
to have equal as at the various redshift Zi, where model spec- 
tra are compared, different models are started with (slightly) 
different ctr at Zin. (Here an is the m.s. fluctuation over the 
scale R; e.g.: as = ""sh-iMpc-) The smallest resolved scale 
Rr yields the greatest relevant aji^ . The selected Zin is such 
that residual non-linear effects cannot induce discrepancies 
among models, over the scale Rr, exceeding 1:10* (Crocce 
et al. 2006). For the ACDM model (DMO and DMG) we 
also ran a smaller box L — 64/i~^Mpc (ACDM-small), to 
examine the resolution effects. 

One of the key analyses in this work are power spectra 
of the matter density field. The matter density field is found 
by interpolating the particle distribution onto a regular Ng x 



Figure 3. Ratio between stellar mass and total mass in haloes. 
Empty and solid squares show the results for the large box (L = 
256) in the ACDM and dDE models. Crosses show results for the 
ACDM small box (L = 64). The solid line represents the recent 
observational results from Moster et al. (2010), dotted lines show 
a 40% scatter (~ 2(t) around the mean. Notice that dDE yields a 
smaller Mstar, at a fixed Mb,aio mass. The slightly overproduction 
of stars at high halo masses is due to the absence of AGN-like 
feedback in our hydrodynamical simulations. 



Ng X Ng grid using the Cloud-in-Cell algorithm. The power 
spectra are calculated using a EFT (Fast Fourier Transform) 
of the matter density field. For our analysis we use Ng = 
2048 for the large boxes (L = 256/i"^ Mpc) and Ng = 512 
for ACDM -small, in order to cut the spectra at the same 
frequency. 

Another key piece of analysis was finding haloes. We 
identified collapsed structures using the spherical overden- 
sity (SO) algorithm. We used a time varying virial density 
contrast from the fitting formulae of Mainini et al. (2003). 
The halo catalogue includes all structures with more than 
200 DM particles. For DMG simulations we estimated the 
halo masses by taking into account all particles inside the 
virial radius (see Maccio et al. 2008 for further details on 
our halo finding algorithm). 



5 MASS FUNCTIONS AND STAR 
FORMATION 

A significant test of our simulations are halo mass functions. 
Figure [2] compares the cumulative mass function A''(> M) 
in our ACDM simulations with the Sheth & Tormen (2002) 
predictions for z = and z = 1. Figure [2] shows the haloes 
found in the DMO and DMG simulations for both 256 and 
64/i~^Mpc boxes. The consistency between the DMO and 
DMG simulations is readily apparent, as is the continuity of 
the mass function over 4 orders of magnitude. 

We make a more direct comparison with observations in 
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Figure 4. Star formation rate in DMG simulation. The upper 
panel shows the evolution of the star production rate vs. z: solid 
histograms are ACDM and dDE rates (blue and red respectively) 
in the 256/i~ Mpc box; the dashed (black) line is Hopkins & 
Beacom fit (2006, assuming a Salpeter IMF); all curves are nor- 
malized to the z = value. The apparent high— 2; discrepancy is 
due to the low-mass halo cutoff (« 200mc ~ 1.3 X 10^-^/i~^Mq) 
in the simulations. The lower panel shows the ratio between star 
productions in dDE and ACDM models. The former model clearly 
exhibits a lower stellar production at all redshifts. 



Figure |3] which shows how the ratio Mstar/Mhaio depends 
on the halo mass in both 64 and 256/i~^Mpc boxes. Figure 
[3] shows the average Mstar/Mhaio for each mass bin. These 
are compared with the predictions from the halo occupation 
model of Moster et al. (2010). Each simulation point is ob- 
tained from equal halo numbers and horizontal error bars 
yield the m.s. spread of \og{Mhaio). Even though the simu- 
lations are run at low resolution relative to star formation 
scales, the number of stars formed in the haloes generally 
follows the halo occupation distribution trend. It is possible 
that too many stars form in haloes in the highest M^aio bin. 
This is likely due to overcooling that plagues hydrodynamic 
simulations and might be fixed by including AGN effects in 
the simulations. 

Most important to the current study, Figure |3] shows a 
significant reduction of Mstar/Mhaio from the ACDM simu- 
lations to those run with dDE. Such reduction is significant 
and one should not be misled by the apparent consistency 
within error bars: they are large because of the spread of 
Mstar/Mhaio through individual haloes, but the reduction 
of Mstar/Mhaio IS uot randomly distributed, exhibiting just 
a slight increasing trend from smaller to greater haloes. We 
shall return on this point below. 

Let us however remind that star formation in hydrody- 
namical simulations strongly depends on numerical (mainly 
mass) resolution (e.g. Mayer et al. 2008 and references 
therein). We therefore expect the result shown in Figure 



[3] to change if the same simulations are run at higher resolu- 
tion (not only because we would resolve lower mass haloes) . 
Naively one could expect to see a higher stellar production 
and, therefore, a significant deviation of simulated curves 
from observed ones. The point of our comparison with real 
data is not to show that our resolution allows intrinsically 
correct results, but that, with resolution adopted m this spe- 
cific work, hydro simulations produce a reasonable amount 
of stars. 



Figure U shows the cosmic star formation history in the 
256/i"^Mpc boxes for both the ACDM and dDE simula- 
tions. The top panel shows the simulated star formation 
rate density (SFR, p*) normalized to the observed z = Q 
value from Hopkins & Beacon (2006). Due to the low res- 
olution of the simulations, they underproduced stars by a 
factor of ~ 1.5. We renormalized these values to show that 
the shape of the star formation history is similar to obser- 
vations. The observational results are shown with constant 
±20 % error bars that roughly approximate the 2-a uncer- 
tainty of their results. Besides the slight deficiency in star 
formation, there are two key differences between the simu- 
lations and observed cosmic SFH. The Hopkins & Beacon 
(2006) observations peak at z ~ 2.3 while the simulations 
peak around 2 ~ 1.5. The onset of star formation is also later 
in the simulations than the observations. This discrepancy 
is again a result of low resolution in the 256/i~^Mpc boxes. 
Christensen et al. (2010) showed that stars don't form in 
halos with less than 200 gas particles. Such haloes have a 
mass of ~ 1.3 x lO^'^ h~^MQ in our simulations. Sufficient 
numbers of these haloes do not form in our simulations until 
2 = 4 to form stars. The SFH of the 64/i~^Mpc box more 
closely follows the observed SFH with a peak around 2 ~ 2.5 
and an early onset of star formation. The results in this case 
are also affected by significant sample variance due to the 
limited size of the simulated volume. 



The bottom frame of Figure [3] shows the ratio between 
p* the dDE and ACDM simulations, highlighting the re- 
duced star formation found using the dDE model. The de- 
crease is ~ 20% and it is worth noting that it originates sim- 
ply from a change of the DE equation of state. The models, 
in fact, have exactly the same age at z = and identical 
cosmological parameters, apart from w. The ~ 20 % dis- 
crepancy shows how different the evolutionary histories and 
concentration distributions are between the two cosmologies. 
The magnitude of the discrepancy is larger than anticipated 
and we will devote future work to determining the cause of 
such a significant difference. 

Before concluding this section let us however return on 
the issue concerning resolution. The comparison with real 
data just shows that our hydro simulations give a reason- 
able star production rate, as function of halo mass, and a 
reasonable distribution of such rate, as a function of redshift. 
This ensures us that the effect of baryonic matter (mainly 
stars) on the total matter power spectrum should be realis- 
tic. This will be true especially on the spatial scales {R) that 
we will address in this work (for fc « 10 — >■ i? « 300/i~^ kpc), 
scales that are much larger than the typical sizes of galaxies. 
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Figure 5. Spectra for different components in the ACDM 
256h~^yipc box. Results for dark matter only (DMO) and dark 
matter & gas (DMG) components are shown using different col- 
ors. 



6 SPECTRA 

Now that we have shown that our hydrodynamic model is 
reasonable on the large scales we wish to study, we turn our 
attention to how the density power spectra depend on the 
DE model used. 



6.1 Hydro vs. N~body simulations 

Figure [5] shows a comparison between DMO and DMG spec- 
tra at 2 = 0. Both the DM and total power spectra of the 
coUisional case show a substantial enhancement especially 
at scales k > 2h Mpc^^. This scale corresponds to the most 
massive structure in our simulation and reflects the fact 
that gas has condensed and cooled inside those structures 
driving a contraction of the haloes themselves (e.g. Gnedin 
et al. 2004). The gas shows a bias with respect to the other 
DMG spectra at large scales k < 8 — lOh Mpc"^. This effect 
has been shown in previous work (Jing et al. 2006, Rudd 
et al. 2008) and is related to gas moving from large scales 
to small scales due to cooling. The cooling is most obvious 
at fc > lO/i Mpc"^, where the gas spectra has more power 
than DMO, reflecting the presence of dense clumps of cold 
gas at small scales, in the center of dark matter halos. 

The differences between DMO and DMG can be better 
appreciated in Figure [S] where the ratio between the different 
spectra is shown as a function of redshift. It is interesting 
to note that the large scale bias of the gas spectrum already 
shown in Figure [S] is strongly redshift dependent and starts 
to develop only after z ~ 1 and is related to gas cooling and 
consumption within dark matter haloes. Another interesting 
feature at z = is the higher power in the DM spectrum 
with respect to the total one in the DMG simulation at 



Figure 6. DMG vs. DMO spectra comparison for the 256 
h~^ Mpc box at different redshifts. 



intermediate scales (2 < k/{h~^ Mpc) < 10). This is due to a 
combination of two effects: the large bias in the gas spectrum 
(which reduces the total one) and the small contribution of 
the stellar component on these scales, that are larger than 
typical galactic scales. For k > lOh Mpc"^ the gap between 
the total and DM spectra vanishes thanks to the increased 
contribution of cold gas and stars. We will come back to 
discuss the effect of stars and resolution on DMG spectra in 
section 16.21 

A further interesting result is the comparison between 
DMO or DMG spectra and the halofit expressions (Smith 
et al. 2003). Only the 256 h~^Mpc box can be used for this 
comparison. The results are shown in Figure [T] The validity 
of HALOFIT up to A; < 2-3/iMpc-^ is confirmed, although 
discrepancies 0(6 %) are found in some spectral ranges. A 
clear conclusion is that, when aiming at 1 % precision, the 
HALOFIT expressions need to be significantly improved. 

6.2 Box size and resolution effects 

Figure |8] compares the power spectra of the large and small 
DMO boxes with the halofit expressions at 2 = 0. While 
the power spectrum from the 25Qh~^Mpc box agrees with 
HALOFIT across the entire range of its validity, the spectrum 
from the 64/i-^Mpc box exhibits a shortage of power in 
the spectral region where non-linearity begins to display its 
effects. This shortage of power is due to the lack of long 
wavelength components in the small box, and clearly shows 
that L ~ 60-70/1"^ Mpc boxes are effected by significant 
sample variance. Thus, we cannot draw quantitative con- 
clusions above k ~ 0.5-1 ft Mpc"'^ from such small boxes. 
Conversely hydrodynamical simulation results could be de- 
pendent on numerical resolution, and the higher resolution 
possible in small boxes make them a useful tool to address 
resolution issues. 
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Figure 7. Comparison between spectra from simulations and 
halofit predictions for the ACDM model in the L = 256/i~^ Mpc 
box. The dashed (black) curve shows the ratio between matter 
spectrum in the DMO simulation and halofit. The solid (red) 
curve shows the same ratio but for the total (dark+star+gas) 
spectrum in the DMG simulation. 
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Figure 8. Comparison between DMO spectra in the large (L = 
256h~^Mpc, solid line) and small (L = 64/1"^ Mpc, dotted- 
dashed line) boxes vs. HALOFIT expressions for the ACDM cos- 
mology. The large box is in reasonable agreement with HALOFIT 
in its validity range (fc < 3feMpc~^); while the small box shows a 
shortage of power in that range, due to the lack of long wavelength 
components. 



Figure |9] shows the fractional difference between DMO 
and DMG simulations in the small box. The gas component 
shows a slightly different behaviour than it did in Figure |6l 
In the small box, the bias at large scales is already present 
at z = 1.2 and grows substantially to 2 = 0, propagating its 
effects at smaller and smaller scales. This is a consequence of 
the enhanced SF in the small box that, as mentioned in ^ 
peaks at z = 2.4. Owing to the (smaller) discrepancies in the 
overall spectra, it seems natural to argue that a systematic 
power depletion in the gas spectrum occurs because of stars 
formation in small scale sites. This is confirmed by a compar- 
ison with Rudd et al. (2008) spectra, shown in their Figure 
2, right-hand panels. In their case, the power depletion in 
the gas spectra is even stronger, so that the overall spectra 
are systematically above CDM. However, Rudd et al. (2008) 
state in their paper that their simulations are characterized 
by an excess in star formation by a factor ~ 3. Because of 
the high star formation, they see the same depletion in their 
gas power spectra, only stronger. The enhanced star forma- 
tion in the small box also has an effect on the ratio between 
total and DM spectrum in the DMG simulation. In the small 
box, the total spectrum is higher than the DM spectrum for 
k > lO/i Mpc~^ thanks to the stellar component. 

The effect more important than resolution is the intrin- 
sic lack of power in the small box. Thus, the total spectrum 
is strongly reduced in the small box with respect to the large 
one. This difference is on the order of 50% at fc = lO/i Mpc~^ 
and even larger, up to 100% for k — 20h Mpc~^ (see also 
Levine & Gnedin 2006). It is interesting to note that our 
results on the total, gas and DM spectra behaviour in the 
large box are closest to the results of Jing et al. 2006 (who 
used a box of 100 h~^ Mpc ), who also found the DM spec- 
trum to lie above the total one for intermediate values of 
k. 

The two conclusions we draw from this comparison be- 
tween the small and large box are: (i) L ~ 60-70 /i~^Mpc 
boxes are effected by significant sample variance, frequently 
yielding potentially misleading results, and hardly allowing 
quantitative conclusions above k ~ 0.5-1 /iMpc"'^; (ii) in 
turn, boxes with a side L ~ 250-260 /i~^ Mpc yield results 
that will only require slight quantitative adjustments (pos- 
sibly at higher resolution) for k > 6-7/iMpc~^. 

We are thus left with no reasons why we should not 
trust spectral comparisons between different models, shown 
in the forthcoming Figures. If some (minor) bias is present, 
it should affect related spectra in a similar fashion. 

6.3 dDE and spectral regularities 

We now turn our attention from the effects of resolution and 
baryons on power spectra to the effects of using different 
cosmologies. 

Figure [TO] shows the fractional spectral differences be- 
tween ACDM and dDE (upper panel) and ACDM and 
ACDM 1 (lower panel). Two main effects are apparent: (i) 
the curves are much smoother for dDE than for ACDM 1. In 
particular, in the top panel, the discrepancy at z = almost 
vanishes, for the total spectrum. This is because ACDM is 
the yV{z = 0) model of dDE, so these findings are expected. 
In particular, the discrepancies up to ~ 3 % at higher z 
affirm the findings of Francis et al. (2007). The discrepan- 
cies between ACDM and ACDM 1 (bottom panel) are far 
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Figure 9. DMG vs. DMO spectra comparison for the small box 
(L = 256/i~^ Mpc) at different redshifts. 



more exaggerated and are typical of models lying ~ 2(7 's 
away from WMAP5 data. This shows that spectral analysis 
that can detect 1 % differences will successfully discriminate 
among models within the 2-0" curves in Figure [T] (ii) The 
ACDMl model re-affirms that DMO-DMG spectral discrep- 
ancies become most significant at fe ~ 2-3 /iMpc"'^. 

In Section 2 we reviewed the technique presented in 
Paper I, that we shall explicitly test here. We start from the 
assigned model A, which is our dDE model, and find the 
auxiliary models y\^{z) for at 2 = 0, 1 and 2 . They must 
share the values of ijJb,c,m = ^b,c.mh? and ag with A. Let 
us recall that the former request is easily fulfilled, as H^ — 
{8nG/3)pcr and, by multiplying both sides of this relation by 
Qrn. (or fib, fie) we have that ujrn oc p,„ « a~^, independent 
of the model. Conversely, the evolution of as depends on 
the state equation of DE and its values at 2 = and 2 = 24 
can only be worked out once we know the constant DE state 
parameters w{z) of the W{z) models for z = 0, 0.5, 1, 2 . 
The final requirement, causing the dependence on 2 of the 
constant ui's is that w is tuned so that W{z) and A have 
equal distances between 2 and the LSB. 

In Figures [TT] and \T2\ we compare spectra at 2 = and 
0.5, 1, 2 . Notice the extreme expansion of the ordinate units 
in these plots. Besides the gas case, the whole ordinate range 
is within ±1 %, for k < lO/i Mpc"^ 

The basic result is that regularities persist when gas 
dynamics is important, although residual discrepancies (of 
the order of a few permils) are greater in DMG simulations. 
The largest discrepancy is found in the gas spectra, even if 
the gas spectrum distortion at 2 = 0, attaining ~ 5%, ap- 
parently does not imply a significant distortion in the global 
spectrum. Altogether it is fair concluding that even the in- 
clusion of hydrodynamics keeps the discrepancies between 
the auxiliary model and the true model in the permil range 
and this is one of of the major results of this work. 
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Figure 10. Fractional spectral differences comparisons: ACDM 
vs. dDE is in the upper plot and ACDM vs. ACDM 1 is in the 
lower plot in DMO and DMG simulations. Four different redshift 
values are shown. Color meaning are the same as in Figure [5] 



6.4 Impact of the results on next Aveak lensing 
data survey 

Matter density spectra P{k, 2) provide a basic link between 
observational data and theoretical models of the nature of 
DE. We envisage comparing the models with data from to- 
mographic shear surveys. Such surveys allow one to work out 
the angular power spectrum P^{£), for weak lensing conver- 
gence between the i-th and j-th tomographic beams cover- 
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Figure 11. Comparison between A and W spectra, for DMO and 
the various components in DMG. Top and bottom panel show 
results for 2 = and z = 0.5 respectively. 



Figure 12. Same as figure [TTI but for z = 1 (top panel) and 2 = 2 
(bottom panel). 



ing the redshift intervals A; and Aj . It is then easy to show 
that Pij{l) is directly related to P{k,z), being 



;« = (f) 



dzP 



r{0,z 



■,z 



W.{z)WAz)^y{^) 



Here H{z) is the Hubble parameter at the redshift z, r{z, z') 
is the comoving angular diameter distance between z and z' , 
while 

W^{z) = 1.5 n,^{l + z) dung{u)r{u,z)/r{0,z) , (5) 

ng(z) being the comoving galaxy number distribution on 
redshift. Using eq. Q, tomographic shear survey data, for 



100 < I < 3000, will yield matter fluctuation spectra for 
0.3 ~< /c//iMpc-^ <~ 10. 

The essential point, here, is that such spectra arise from 
any possible gravitational source. Any link with assumptions 
on the mass/luminosity ratio is therefore broken. Further- 
more, the expected accuracy of lensing surveys will allow to 
recover P{k,z) with a precision 0{1%). 

For gravitational lensing shear predictions, at first or- 
der, there seems to be no need to take into account the 
baryon component. Baryons are comprise a small fraction 
of the Universal mass fraction compared to DM, and their 
distribution can be only slightly different from CDM. But 
what makes baryons extremely important is the precision 
that shear measurements will reach {0{1%)) in the near 
future. This presents a new challenge for hydrodynamical 
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Figure 13. Evolution of the state parameter wi^z). The A model 
is shown by the solid line and compared with the constant-io val- 
ues of the auxiliary models W (dotted line) computed at different 
redshifts. Solid squares show the w values used in this work. 



simulations. Formerly, a slight inaccuracy in baryon physics 
was a second order problem compared with the high preci- 
sion achieved in tracing the total matter distribution. While 
baryons are less massive than CDM, the way CDM reacts 
to the evolution of the baryon distribution makes baryons 
an essential ingredient. 

Let us then outline a very crucial issue concerning the 
extrapolation of w^z) from tomographic cosmic shear data. 
This point, already made in Paper I, is strengthened by this 
work, taking into account baryon physics. The point is that 
the conversion between the observed evolution of w(2:) and 
the true one of the underlying model is not a straightforward 
operation. 

Figure [13] shows the redshift evolution of the state pa- 
rameter of the A model (solid line) compared with the values 
of the constant-w in the different auxiliary models (W) at 
different redshifts (solid squares). The figure outlines that 
cosmic shear measures, fitted assuming constant w in each 
z-hin (i.e. the solid squares) would infer an evolution of w 
with redshift (the dotted line) radically different from the 
true underlying model A (solid line). 

The mapping between the w values best fitting data 
in the tomographic layer about z and the real evolution of 
w{z) should be done with extreme care as shown in Casariru 



(2010). 



7 DISCUSSION AND CONCLUSIONS 

In this work, we employed large scale hydrodynamical sim- 
ulations, including gas cooling, star formation and Super- 
novae feedback. Our simulations perform quite well when 
compared to observational results. They produce a sizable 
stellar component in collapsed dark matter haloes and star 
formation histories comparable with observations. We are 



aware that these results strongly depend on the adopted 
(low) resolution, nevertheless the similarity to observations 
makes us confident that the effect of baryons on the mat- 
ter power spectrum in our simulation is sufficiently real- 
istic. That said, further work is needed to better address 
the impact of the physics of star formation on the final re- 
sults and, most important, to provide better tests of sam- 
ple variance. In particular, this work makes clear that no 
really quantitative predictions, unaffected by sample vari- 
ance, can be drawn from simulations, in i ~ 50fe~^ Mpc 
boxes. Matching observed accuracies of C(l%) will require 
boxes a few hundred Mpcs on a side. In principle a box 
with L = 256/i~^Mpc, as the one we mostly used here, 
should be adequate to delve into haloes, providing spectra 
up to fc ~ 10/iMpc~^, with the required accuracy. However, 
the spectra obtained from our large box differed in a non 
negligible way from the ones obtained from a smaller box 
(L = 64/i~^ Mpc). Although being fully compatible with 
the sample variance of the smaller box, this suggests further 
investigation will be required. 

These residual uncertainties do not prevent us, though, 
from testing the efficiency of a method introduced in Paper 
I, that allows us to work out power spectra prediction for any 
DE model based on models with a constant equation of state 
iw). We found up to z = 2, that the auxiliary W models 
(with constant w) yield spectra consistent with our assigned 
A model (with a polynomial Vli{z)) within a few permils. 
The test was performed for a single A dDE model. We plan 
to extend such a test to more dDE models in forthcoming 
work. There is however no reason to expect that the point 
selected on the Wo-w' plane owns peculiar properties. 

The residual discrepancies between A and W models 
exhibit a specific trend, already visible in Paper I: discrep- 
ancies (of a few permils) are greater for smaller z values. 
Also DMG results follow the same rule, which is even even 
more pronounced for baryon spectra. Gas spectra discrep- 
ancies often exceed the permils but are systematically com- 
pensated by other components. The success of our technique 
was not unexpected, but discrepancies at such a tiny level, 
also for k > 2-3/iMpc~^, exceeded our expectations. 

The requirement that two models have the same ujc.b 
and crs{z) introduces basic similarity: they will have similar 
baryon acoustic oscillations and typical fluctuation ampli- 
tudes. Then, the requirement that w yield the same confor- 
mal time r elapsed since recombination for assigned and 
auxiliary models apparently sets an equal timing for the 
gradual onset of the non- linear evolution on analogous struc- 
tures. However, while r is a fair time coordinate for struc- 
tures slightly detached from the cosmic background, it no 
longer plays such a role in cosmic islands, which have aban- 
doned the cosmic expansion and evolve on their own, for- 
getting the continuing increase of the scale factor a . Within 
these structures, space-time is substantially Minkowskian, 
and time evolution is set by the ordinary time t . 

These are, however, two extreme situations. The scales 
which "recently" abandoned the overall expansion and en- 
tered a non-linear regime are in an intermediate position. 
Our inspection mostly concerns such scales. Accordingly, the 
success of an equal-r requirement is anything but granted. 

Here we use hydrodynamical simulations to study a 
larger k domain (compared to Paper I), where the effects of 
baryons cannot be neglected anymore. By introducing non- 
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gravitational efTects we expect dynamics to be regulated by 
the ordinary time t, thus producing larger discrepancies at 
scales k > 2-3 /iMpc~^. Figure [TT1 rejects this hypothesis 
by showing that the matter distribution is still ruled by the 
primeval clock. In particular, for k > 2-3/iMpc~^, the over- 
all DMG spectral discrepancies, at z = 0, are perhaps even 
smaller than DMO. Curiously enough, the opposite seems 
true at some greater z's, suggesting that most discrepan- 
cies found are pure noise, rather than systematics. Clearly 
for larger fc's, when we delve into regions that virialized long 
ago, A and W models unavoidably differ. However, while we 
restrict ourselves to the scales that forthcoming weak lens- 
ing surveys will study, the Paper I requirement still works 
unexpectedly well. 

One ought to beware, however, that the success of our 
approach does not imply that constant-^; and variable-w 
models are indiscernible. What makes them observationally 
different is the measured dependence of w on the redshift z. 

An essential issue, however, is that the observed constant 
state parameter w{z) will NOT be the physical one of the DE 
component. Although we expect data to fit, at any z, the 
same density and Hubble parameter values, the observed 
'w{z) will be the state parameter that a constant- w model 
must have, at such z, to grant the same distance between 
there and the LSB. The relation between the observed w 
and the intrinsic one, for the DE component, should be put 
under control; otherwise there is a severe risk that future 
data analysis lead to DE state equation misinterpretation. 

A final point we wish to make is that halofit ex- 
pressions, so useful for predicting spectra until now, have 
been tested and found compatible with their claimed accu- 
racy. Unfortunately, accuracy at (~ 6 %) no longer matches 
our needs. Furthermore, among models that approach the 
present data, halofit expressions only need to handle 
ACDM . If one aims to improve halofit expressions, there 
is no longer any need to cover the wide array of models that 
once were relevant, but now are far from observation. 

Amelioration, however, is needed in two separate direc- 
tions: (i) The spectral range covered by predictions must 
reach k ~ 10/iMpc~^, which requires that gas be consid- 
ered; Rudd et al. (2008) had actually attempted an exten- 
sion in this direction, but limited themselves to ACDM cos- 
mologies, (ii) DE state parameters different from constant 
w = — 1 should be included. Moreover, they should aim for a 
precision of Oil %), which presents a severe challenge, even 
for a much more restricted model range. Nevertheless when 
HALOFIT expressions of the required precision will be avail- 
able, even only for constant-w models, thanks to the results 
of this work it will be straightforward to extend them to any 
w{z) law. 



ACKNOWLEDGEMENTS 

Giuseppe Murante is thanked for fruitful discussions. The 
authors acknowledge J. Wadsley for development of the 
GASOLINE code and thank him for its use in this work. Nu- 
merical simulations and spectral analysis were performed 
on the PIA and on PanStarrs2 clusters of the Ma:x-Planck- 
Institut fiir Astronomic at the Rechenzentrum in Garching 
and at the CINECA-Bologna computer center, within the 
frame of the INAF-CINECA agreement (2008-2010). The 



support of AST (Italian Space Agency) through the contract 
1/016/07/0 COFIS is acknowledged. 



REFERENCES 

Brook C. B., Kawata D., Gibson B. K., Flynn C, 2004, 

MNRAS, 349, 52 
Casarini L., Maccio A. V., Bonometto S. A., 2009, JCAP, 

3, 14 (Paper I) 
Casarini L., 2010, NewA, 15, 575 

Christensen C., Quinn T., Stinson G., Bellovary J., Wads- 
ley J., 2010, arXiv, arXiv:1005.1681 
Crocce M., Pueblas S., Scoccimarro R., 2006, MNRAS, 373, 

369 
Dodelson S., Zhang P., 2005, PhRvD, 72, 083001 
Dolag K., Bartelmann M., Perrotta F., Baccigalupi C, 

Moscardini L., Meneghetti M., Tormen G., 2004, A&A, 

416, 853 
Francis M. J., Lewis G. F., Linder E. V., 2007, MNRAS, 

380, 1079 
Gnedin O. Y., Kravtsov A. V., Klypin A. A., Nagai D., 

2004, ApJ, 616, 16 
Governato F., et al., 2004, ApJ, 607, 688 
Governato F., Willman B., Mayer L., Brooks A., Stinson 

G., Valenzuela O., Wadsley J., Quinn T., 2007, MNRAS, 

374, 1479 
Guillet T., Teyssier R., Colombi S., 2010, MNRAS, 405, 

525 
Haardt F., Madau P., 1996, ApJ, 461, 20 
Hearin A. P., Zentner A. R., 2009, JCAP, 4, 32 
Hopkins A. M., Beacom J. F., 2006, ApJ, 651, 142 
Hu W., Tegmark M., 1999, ApJ, 514, L65 
Huterer D., Takada M., 2005, APh, 23, 369 
Jing Y. P., Zhang P., Lin W. P., Gao L., Springel V., 2006, 

ApJ, 640, L119 
Katz N., 1992, ApJ, 391, 502 
Kennicutt R. C, 1998, ApJ, 498, 541 
Klypin A., Maccio A. V., Mainini R., Bonometto S. A., 

2003, ApJ, 599, 31 
Komatsu E., et al., 2009, ApJS, 180, 330 
Kristiansen J. R., La Vacca G., Colombo L. P. L., Mainini 

R., Bonometto S. A., 2009, arXiv,|irXw^:0902.2737 
La Vacca G., Kristiansen J. R., Colombo L. P. L., Mainini 

R., Bonometto S. A., 2009, JCAP, 4, 7 
La Vacca G., Bonometto S. A., Colombo L. P. L., 2009, 

NewA, 14, 435 
Levine R., Gnedin N. Y., 2006, ApJ, 649, L57 
Linder E. V., Jenkins A., 2003, MNRAS, 346, 573 
Maccio A. V., Quercellini C, Mainini R., Amendola L., 

Bonometto S. A., 2004, PhRvD, 69, 123516 
Maccio A. V., Dutton A. A., van den Bosch F. C, 2008, 

MNRAS, 391, 1940 
Mainini R., Maccio A. V., Bonometto S. A., Klypin A., 

2003, ApJ, 599, 24 
Manera M., Mota D. F., 2006, MNRAS, 371, 1373 
Mayer L., Governato F., Kaufmann T., 2008, ASL, 1, 7 
Moster B. P., Somerville R. S., Maulbetsch C, van den 

Bosch F. C, Maccio A. V., Naab T., Oser L., 2010, ApJ, 

710, 903 
Mota D. F., 2008, JCAP, 9, 6 



© 2002 RAS, MNRAS 000, [T]-? 



12 Casarini et al. 

Refregier A., Artiara A., Kitching T. D., Rassat A., 

Scaramella R., Weller J., Euclid Imaging Consortium f. t., 

2010, arXiv, arXivilOOLOMT 
Rudd D. H., Zentner A. R., Kravtsov A. V., 2008, ApJ, 

672, 19 
Sheth R. K., Tormen G., 2002, MNRAS, 329, 61 
Smith R. E., et al, 2003, MNRAS, 341, 1311 
Solevi P., Mainini R., Bonometto S. A., Maccio A. V., 

Klypin A., Gottlober S., 2006, MNRAS, 366, 1346 
Stadel J. G., 2001, PhDT, 
Stinson G., Seth A., Katz N., Wadsley J., Governato P., 

Quinn T., 2006, MNRAS, 373, 1074 
Wadsley J. W., Stadel J., Quinn T., 2004, NewA, 9, 137 
Wetterich C., 1988, NuPhB, 302, 668 
Zentner A. R., Rudd D. H., Hu W., 2008, PhRvD, 77, 

043507 



© 2002 RAS, MNRAS 000, [l]-?? 



